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Strong-field physics is currently experiencing a shift towards the use of mid-IR driving wavelengths. 

This is because they permit conducting experiments unambiguously in the quasi-static regime and 
enable exploiting the effects related to ponderomotive scaling of electron recollisions. Initial mea¬ 
surements taken in the mid-IR immediately led to a deeper understanding of photo-ionization and 
allowed a discrimination amongst different theoretical models. Ponderomotive scaling of rescatter¬ 
ing has enabled new avenues towards time resolved probing of molecular structure. Essential for 
this paradigm shift was the convergence of two experimental tools: 1) intense mid-IR. sources that 
can create high energy photons and electrons while operating within the quasi-static regime, and 2) 
detection systems that can detect the generated high energy particles and image the entire momen¬ 
tum space of the interaction in full coincidence. Here we present a unique combination of these two 
essential ingredients, namely a 160 kHz mid-IR source and a reaction microscope detection system, 
to present an experimental methodology that provides an unprecedented three-dimensional view of 
strong-field interactions. The system is capable of generating and detecting electron energies that 
span a six order of magnitude dynamic range. We demonstrate the versatility of the system by 
investigating electron recollisions, the core process that drives strong-field phenomena, at both low 
(meV) and high (hundreds of eV) energies. The low energy region is used to investigate recently 
discovered low-energy structures, while the high energy electrons are used to probe atomic struc¬ 
ture via laser-induced electron diffraction. Moreover we present, for the first time, the correlated 
momentum distribution of electrons from non-sequential double-ionization driven by mid-IR pulses. 

PACS numbers: 32.80.-t, 32.80.Fb, 34.80.Bm 


Strong-field physics (SFP) is concerned with the ef¬ 
fects related to the interaction of intense electric fields 
with matter. Here, ‘intense’ means that the electric field 
strengths are non-negligible compared to the binding 
fields within matter. Since these field strengths are most 
easily attained with intense ultrashort laser pulses, most 
of the effects and observations stem from the interaction 
with oscillating laser fields. In parallel with the advance¬ 
ment of ultrafast laser science, SFP has developed into 
a mature field of research that is now capable of track¬ 
ing electronic and structural dynamics on the atto- to 
few femtosecond timescales HHZ]. This advance of SFP 
and attoscience has, in turn, generated an upsurge in 
the development of ultrafast mid-IR laser sources due to 
the possibilities when driving strong-field recollision with 
long wavelengths mm- Mid-IR sources present many 
benefits compared to the ubiquitous 0.8 /jm radiation of 
Ti:Sapphire based laser systems. Some of the most rel¬ 
evant to this study are: 1) the unambiguous creation 
of interaction conditions that are conducive for classi¬ 
cal interpretations of experimental results; 2) strong-field 
recollision at low peak intensity in order to avoid ap¬ 
preciable ground state depletion; and 3) the ability to 
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create high energy recollision electrons for imaging ap¬ 
plications. However, long wavelength (ponderomotive) 
scaling comes at the cost of a dramatically reduced sig¬ 
nal El EMU, which translates into reduced statistics or 
long data acquisition times. Exacerbating this problem 
is the fact that traditional electron detection techniques 
only detect a small fraction of the entire momentum 
space (see Sec. IV). These points are often overlooked 


but they present a significant roadblock to the advance¬ 
ment of SFP. This is because experiments will eventually 
become untenable due to the required acquisition times 
exceeding the stability of the apparatus. Here we present 
a methodology that overcomes these limitations by com¬ 
bining two tools that are perfectly suited for mid-IR SFP: 
1) an intense and high repetition rate optical paramet¬ 
ric chirped pulse amplification (OPCPA) based mid-IR 
(3.1 /zm) source [IS] that creates ionization conditions 
deep in the quasi-static (tunneling) regime and that can 
generate electron energies above 1 keV; and 2) a reaction 
microscope (ReMi) detection system that can image the 
three-dimensional (3D) momentum space of the interac¬ 
tion in full particle coincidence. The apparatus is capable 
of generating and detecting electron energies that span 
a six order of magnitude range and thus provides an un¬ 
precedented 3D view of strong-field interactions. 

The paper is structured as follows: We first discuss 
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relevant aspects of strong-field ionization (SFI) (Sec. [I]) 
when scaling the driving wavelength of the radiation and 
not the peak intensity; the metric we use to describe 
such scaling is the adiabaticity, or ‘Keldysh’, parameter 
(Sec. 0 . Next, the validity of the dipole approximation is 
discussed in the context of wavelength scaling (Sec. III). 
We then highlight the critical role of detection systems 
with specific attention paid to molecular targets and the 
importance of measuring all ionization fragments in full 
coincidence (Sec. IV). We proceed to demonstrating the 
versatility of ou r m ethodology by investigating electron 
recollision (Sec. VI) at the extremes of a 550 eV wide 


electron spectrum. We highlight two cases: In Sec. VIA 


the low-energy part of the spectrum (meV) reveals the 
recently discovered low energy features for both single 
and double ions while in Sec. VIB the high-energy por¬ 
tion (hundreds of eV) is used to extract accurate differen¬ 
tial cross-sections (DCSs) over a wide range of recollision 
energies using laser-induced electron diffraction (LIED). 
Lastly, we discuss our methodology in the context of the 
current status and expected direction of SFP research. 


I. STRONG-FIELD IONIZATION 
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FIG. 1: Scaling of the Keldysh parameter versus driving 
wavelength for xenon atoms ionized with an intensity of 
1 x 10 14 Wcm -2 . The regions of 7 1 and 7 <C 1 are 

highlighted by grey and red shading, respectively. The 
ubiquitous Ti:Sapphire wavelength of 0.8 fim (red 
square) and our OPCPA outputs (red circles) are 
presented. The inset shows the kinetic energy E kin of 
elastically back-scattered electrons (9 T = 180°) versus 
their return energy E T (for long and short trajectories) 
according to the semi-classical ‘recollision model’ in 
units of the ponderomotive energy Up. 


Strong-field interactions can be qualitatively under¬ 
stood by the semi-classical recollision model pj5|, which 
describes the motion of an electron after strong-field in¬ 
duced tunnel ionization. According to this model the 
electron is emitted close to the peak of the oscillating 


laser electric field, after which it accelerates within the 
field before re-colliding with its parent ion roughly three- 
quarters of an optical cycle later. In the recollision model 
returning electron trajectories can be distinguished as 
long and short depending on the time they spent ac¬ 
celerating in the held before returning to their parent 
ion. The maximum return energy of E r max = 3.17 Up 
(Up being the ponderomotive energy experienced by the 
detached electron) is achieved when the electron is emit¬ 
ted at a phase of about 17° after the held maximum. 
This results in the electron returning to the ion at about 
255°. Electrons emitted before 17° return after 255° and 
are considered long trajectories, while those emitted after 
17° return before 255° and are considered short trajec¬ 
tories. The model typically 1) neglects depletion of the 
ground state and Stark shifts; 2) assumes that only one 
electron participates; and 3) omits the influence of the 
Coulomb held on the electron’s trajectory. Even with 
such assumptions, numerous potential outcomes of the 
electron-ion re-collision can be qualitatively described: 
the returning electron may re-combine with the ion and 
emit a high energy photon in a coherent frequency up- 
conversion process known as high harmonic generation 
(HHG); it may inelastically scatter of the ion and cause 
non-sequential double (or multiple) ionization (NSDI); or 
it may elastically scatter off the ion resulting in a hnal 
kinetic energy E k i n that is determined by the energy E T 
at which it returns to the parent ion and the angle 8 r 
at which it scatters. The inset of Fig. [I] shows the de¬ 
pendence of Eirin on E r at d r = 180° for long and short 
electron trajectories. A framework to model the ioniza¬ 
tion step of such strong-field interactions was proposed 
in Keldysh’s seminal work from 1965 [20]. The dimen¬ 
sionless adiabaticity parameter 7 that is commonly used 
to describe the interaction conditions was introduced, 


7 \l~2up~ Vix' (1) 

The ‘Keldysh’ parameter 7 is a function of the pondero¬ 
motive energy Up oc IX 2 , where I is the laser peak inten¬ 
sity and A is the central driving wavelength. The scaling 
of the Keldysh parameter as a function of the driving 
wavelength is presented in Fig. [T] for a typical peak in¬ 
tensity of 1 x 10 14 W cm -2 interacting with a xenon (Xe) 
atom (Ip ss 12 eV). A rapid decrease of the Keldysh 
parameter can clearly be observed as the driving wave¬ 
length is increased from the visible (A < 0.8 /im) through 
the near-IR (0.8 /cm < A < 3/rrn) and into the mid-IR 
(A > 3/cm) [SI- The circular data points represent the 
wavelengths that are generated by our OPCPA source. 
The conventional interpretation of 7 is as a metric to 
distinguish between the extremes of two different descrip¬ 
tions of ionization: 1) the so-called ‘multi-photon’ (MP) 
regime where 7 1 and experimental electron kinetic 

energy distributions generally feature peaked behavior; 
and 2 ) the so-called ‘tunneling’, or to be precise ‘quasi¬ 
static’ (QS), regime where 7 <C 1 and the distributions 


















3 


typically exhibit continuous behavior. The case when 
7 sa 1 is known as the ‘transition regime’ and ionization 
is characterized by features from both MP and QS ion¬ 
ization 1221 . The square data point in Fig. [I] is within 
this region and it serves as a reference for the ubiquitous 
Tirsapphire radiation at 0.8 /im. Most SFI experiments 
performed in the past couple of decades have utilized 
laser systems based on this technology, therefore theo¬ 
retical understanding has evolved in this context. 

Several motivating factors exist for investigating SFP 
in the QS regime. The most obvious is the desire to 
match the assumptions of theoretical models to further 
our understanding of SFP in general. Another important 
aspect is the fact that QS conditions are a prerequisite 
for invoking classical recollision models, which are the 
cornerstone for developing attosecond science or recolli¬ 
sion based time resolved imaging. The requirement for 
unambiguous conditions was strikingly demonstrated by 
ionizing atoms with long wavelength radiation that re¬ 
sulted in unexpected structures at low electron kinetic en¬ 
ergies [SII23]- These structures were not predicted by the 
most widely used approximate theoretical models includ¬ 
ing the workhorse of the field, the strong-field approxi¬ 
mation (SFA). As these models are utilized to interpret 
many experimental results, including attosecond science 
experiments, it was worrying that such obvious features 
could evade prediction and detection for so long. It was 
not until improvements were made to both the SFA and 
to classical models that multiple electron re-scattering 
due to an interplay of the laser and Coulomb fields could 
be identified as the responsible mechanism . Since 

these initial results a number of other lower energy fea¬ 
tures have also been observed for both atoms and small 
molecules during long wavelength investigations within 
the QS regime HU -'T3 2U and the responsible mecha¬ 
nisms have subsequently been identified |35l . Moreover, 
the description of the unbound electron being dominated 
by the laser field within the QS regime is a prerequi¬ 
site to interrogate inelastic re-scattering processes. The 
mechanisms behind multiple NSDI connected to electron- 
electron correlations have been extensively studied and 
theoretically modeled for near-infrared and high inten¬ 
sity fields, yet questions remain whether the findings hold 
for different regimes [36] • Finally, the LIED technique, 
which has recently been utilized to image the structure 
of both diatomic m and polyatomic ED molecules, re¬ 
quires invoking the classical recollision model. Here the 
influence of the strong field is removed and the return¬ 
ing electron wavepacket is treated as a plane wave upon 
recollision. These motivating factors make it clear that 
performing SFI experiments in the QS regime is advan¬ 
tageous in many ways. 


II. SCALING OF THE KELDYSH PARAMETER 

The most straight forward method to reduce 7 is to 
increase the laser peak intensity I. Unfortunately, from 



FIG. 2: The ionization fractions calculated at the peak 
of a six cycle pulse for targets with different ionization 
potentials while keeping the Keldysh parameter 
constant at 7 = 0.3. The squares represent 0.8 /im 
radiation while the circles are the wavelengths of our 
OPCPA. 


an experimental perspective, this approach is severely 
constrained by the ionization potential Ip of the target 
because the ground state population is depleted when 
the so-called saturation intensity /sat is reached for a 
given pulse duration. Such intensities cause the frac¬ 
tional amount of re-colliding electrons to reduce due to 
increased direct ionization on the leading edge of the 
pulse [29:. This is detrimental for re-scattering based 
experiments, therefore it is beneficial for the ionization 
fraction Ip to be kept much lower than unity. Strictly 
speaking, as /sat sets the highest effective intensity that 
a target can experience, it also sets the lower bound on 
the achievable 7 [39]. In Fig. [ 2 ] the ionization fraction at 
the peak of a six cycle FWHM pulse is illustrated j3D] for 
a target with Ip = 12 eV (dashed blue curve) as a func¬ 
tion of the laser wavelength while keeping the Keldysh 
parameter constant at 7 = 0.3. Such an ionization po¬ 
tential could represent targets such as Xe or acetylene 
(C 2 H 2 ), which have ionization potentials of 12.1 eV and 
11.4 eV, respectively. For 0.8 /mr radiation, a peak in¬ 
tensity of 1.1 x 10 15 W cm -2 is required to reach 7 = 0.3, 
which results in complete ionization of the target before 
the peak of the pulse. Clearly such experimental con¬ 
ditions cannot be used to unambiguously investigate re¬ 
scattering phenomena. The same calculation for 3.1 /mr 
radiation, on the other hand, results in a peak intensity 
of 7.5 x 10 13 W cm -2 and an ionization fraction of only 
4 x 10~ 3 . In fact, it is found that in order to have an ion¬ 
ization fraction < 0.1 at 7 = 0.3, a driving wavelength 
longer than 2.2 /(m is required. This analysis highlights 
the problem of intensity scaling: it is difficult to generate 
QS conditions using near-IR radiation without reaching 
ionization saturation. 

Ionization saturation becomes even more problematic 
when larger polyatomic molecular targets are investi- 
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gated, as they typically have lower ionization potentials 
and therefore ionize more easily. To demonstrate the 
consequences of this aspect, the above calculation was 
repeated for an ionization potential of about Ip = 8 eV 
(solid red curve in Fig. [2]), which could represent a poly¬ 
atomic molecule such as naphthalene (CioHg) [41] . In 
this case saturation persists until a wavelength of about 
2 /im. At 3.1 fj, m, however, the calculated ionization 
fraction is still only 0.16, which is adequate for the in¬ 
vestigation of re-scattering phenomena. Clearly, a more 
preferable and robust method to achieve QS conditions 
is to scale the wavelength rather than the intensity [42j. 
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FIG. 3: The quadratic scaling of the maximum return 
energy with wavelength for an intensity of 
1.0 x 10 14 Wcm -2 . The red square represents 0.8 /im 
radiation while the red circles are the wavelengths of 
our OPCPA. The inset shows an electron energy 
distribution, which extends up to 550 eV, that was 
measured after ionization of C 2 H 2 using the 3.1 /mi 
output. 

The dependence of 7 with A exhibits a favourable 
quadratic scaling of the ponderomotive energy Up. As 
predicted by the recollision model, the time that the elec¬ 
tron accelerates in the field is proportional to the laser 
field period. For this reason mid-IR radiation, which has 
a longer optical period, results in much higher energy 
returning electrons. Increasing the driving wavelength 
from 0.8 /mi to 3.1 /tm, for example, results in a fac¬ 
tor of 15 increase in the ponderomotive energy Up. The 
recollision model also predicts that the maximum energy 
that an electron can have upon return to its parent ion is 
equal to F r ,max = 3.17 Up. The trend of this maximum 
return energy versus the driving wavelength is presented 
in Fig. [ 3 ] for a typical peak intensity of 1 x 10 14 W cm 2 . 
The square data point shows the maximum energy gen¬ 
erated by 0.8 /tm radiation as a reference and the circu¬ 
lar data points represent the return energies obtainable 
with the wavelengths generated by our OPCPA source. 
Return energies of hundreds of eV can be created by scal¬ 
ing strong-field interaction to mid-IR wavelengths, com¬ 
pared to the tens of eV by standard 0.8 /tm radiation. 


These much higher energy electrons can be used as a 
probe of nuclear, rather than electronic, structure and 
dynamics as is discussed in Sec. |VIB| The inset of Fig. [3] 
shows an experimentally measured electron energy dis¬ 
tribution that ranges from 0 — 550 eV. In this case the 
intensity was such that the maximum return energy was 
Fr,ma.x = 175 eV. Note that the recollision model predicts 
elastic rescattering to result in overall electron energies 
of £?kin 10 Up (see inset of Fig. [I]). In Sec. [vT]we show 
that, with the appropriate detection system and analysis 
techniques, the entire bandwidth of such electron spectra 
can be utilized to investigate meV features and to extract 
elastic DCSs at energies of hundreds of eV. High return 
energies also give access to studying inelastic rescatter¬ 
ing processes, such as NSDI, not only from valence elec¬ 
trons but also from inner shells [43] . One may expect 
interesting new insights into the multi-electron correla¬ 
tions in many body systems |36l with mid-IR sources. In 
Sec. I VI A~2l we show first results of mid-IR driven double 
ionization with full coincidence detection. 


III. VALIDITY OF THE DIPOLE 
APPROXIMATION 



FIG. 4: The validity of the dipole approximation as a 
function of the driving wavelength and peak intensity. 
The limit where the magnetic field causes the electron 
motion perpendicular to the laser polarization to 
become greater than a Bohr radius is given by the solid 
black line (/3q = 1.0 a.u.). Also shown are the limits for 
/Jo = 0.5 a.u. (dot-dashed line) and /Jo = 2.0 a.u. 
(dashed line). Most 0.8 /im experiments are performed 
in the ‘valid’ region (red squares) while many near-IR 
and mid-IR experiments are close to the limit (Ref. [TT] 
- purple diamonds, Ref. [46; - orange upwards triangle, 
Ref. 0 - green downwards triangle and Ref. m - white 
cross). The range of intensities available to the two 
outputs of our OPCPA, and therefore the regions over 
which we can probe the dipole approximation, are 
presented by blue circles. 
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Many theoretical models used in SFP utilize the dipole 
approximation |44j . in which the influence of the mag¬ 
netic field component on the electron motion is assumed 
to be negligible compared to that of the electric field 
component. However, when the electron displacement 
perpendicular to the field polarization direction (/?o) be¬ 
comes comparable to a Bohr radius (i.e. /3q ~ 1 a.u.), 
which can occur at long wavelengths and high intensi¬ 
ties [Ml 35] > the validity of the dipole approximation be¬ 
comes questionable. The limit of /3q = 1 a.u. is presented 
in Fig. [4] (solid black line) for a wide range of experimen¬ 
tal driving wavelengths and peak intensities. The regions 
where the dipole approximation is valid (invalid) are rep¬ 
resented by white (red) shading. Also shown are the lim¬ 
its for when /3 0 = 2 a.u. and /? 0 = 0.5 a.u. (dashed and 
dot-dashed black lines, respectively) to highlight the fact 
that the validity of the approximation is not defined by 
a clear threshold. 

Most non-relativistic experiments at the standard 
0.8 /tm wavelength are performed at intensities between 
5 x 10 13 — 5 x 10 15 W cm -2 (red squares and dotted line) 
where the dipole approximation is generally assumed to 
be valid. At longer wavelengths the limit of (3q = 1 a.u. is 
reached at lower intensities and a number of experimen¬ 
tal measurements with conditions near or above this limit 
have already been reported in the literature (Ref. [8] - 
green down-triangle, Ref m - purple diamonds, Ref. [46] 
- orange upwards triangle and Ref. [47] - white cross). 
Also shown are the intensity ranges over which the 1.7 /im 
and 3.1 /im of our OPCPA can be utilized for SFI exper¬ 
iments. The large intensity range of our 3.1 /im radia¬ 
tion opens up the possibility for investigations across the 
f3o < 1 a.u. —► 0o > 1 a.u. transition. 


IV. ON THE IMPORTANCE OF THE 
DETECTION SYSTEM 

The interaction between a strong laser field and mat¬ 
ter is inherently three dimensional (3D). To help visualize 
this scenario, a red ellipse is used to approximate a ficti¬ 
tious 3D momentum distribution in Fig. [5ji. In the case 
of a linearly polarized field, such as considered here, the 
momentum distribution is extended more in the direc¬ 
tion parallel to the laser polarization (longitudinal direc¬ 
tion) than in the transverse directions. An established 
and still popular method to detect such a distribution 
is to use a time-of-flight (TOF) spectrometer where the 
energy of both positively and negatively charged parti¬ 
cles can be calculated based on their arrival time after 
ionization. A typical TOF spectrometer could have an 
acceptance angle of about 0 a = 5° (represented by the 
blue cone), which corresponds to a solid angle that is 
about 10% of the full 47 t sr. For such a detection angle, 
any particle that has a momentum component perpen¬ 
dicular to the TOF axis that is > 9% of its parallel 
momentum component (i.e. p± > 0.09p||) will not be 
detected without rotating the TOF in azimuthal and po- 



(b) Time-of-flight (TOF) spectrometer 




FIG. 5: a) Only certain parts of the full 3D momentum 
distribution (red ellipse) can be measured by a TOF 
spectrometer, however, the laser polarization can be 
varied. After ionization of Xe using intense 3.1 /tm 
radiation, a TOF with an acceptance angle 0 a = 5° 
would measure only a small portion (b) of the full 3D 
distribution collected by the RcMi (c). The highest 
measured electron momentum for this data was between 
6.0-6.4 a.u., which corresponds to energies between 
500-550 eV. 


lar angles about the interaction point. Additionally, any 
angular information within the acceptance cone is aver¬ 
aged over. Small pinholes can be installed to decrease 
the angular resolution but this is obviously accompanied 
by a dramatic drop in count rate. The laser polarization 
can be continuously adjusted between parallel and per¬ 
pendicular (as presented in Fig. [f^i) to the detector axis 
so that different parts of the momentum cloud can be ob¬ 
served. If this rotation is performed numerous times, a 
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2D ‘conical cut’ through the momentum distribution can 
be obtained nn, however, the integration times required 
using this method are long. Recent investigations in the 
QS regime that detected low electron energy (< 1 eV) 
structures where p± ~ P|| [HI I55H55] have highlighted 
the limitations of the TOF measurement technique as 
without the aid of high resolution 3D detection, these 
structures evaded observation. 


An example of the portion of a 3D momentum distri¬ 
bution that would be measured by a TOF spectrometer 
is presented in Fig. §o. This data was produced by tak¬ 
ing the actual 3D momentum distribution, which was 
measured by our ReMi after the ionization of Xe gas us¬ 
ing intense 3.1 pm radiation (Fig. §0. and only keeping 
those counts that were found within the 9\ = 5° cone 
around the laser polarization direction. Only the regions 
of the distributions with p y < 0 are shown in order to 
highlight the higher counts near the origin. The projec¬ 
tions of the three spatial dimensions are also presented 
for both measurement methods. The color bar is valid 
for the 3D cloud but not for the projections as they have 
each been individually rescaled. Clearly, a TOF spec¬ 
trometer only collects a small portion of the full electron 
momentum distribution. Specifically, and as mentioned 
above, any electrons that have a large transverse compo¬ 
nent would not be detected. In fact, for an acceptance 
angle of 9\ = 5° only about 2 % of all electrons are col¬ 
lected. The lack of signal in the TOF projections further 
highlights these issues. We note that the ReMi and TOF 
spectrometer distributions should be compared in a qual¬ 
itative manner only as it is possible for TOF spectrome¬ 
ters to have a different detection angle, to be utilized in 
a stereo configuration, and to be used for varying laser 
polarization directions. 


In addition to momentum distributions being intrinsi¬ 
cally three dimensional, electrons may be detected that 
originate from all atomic and molecular fragments in the 
interaction region. This can include ions and isotopes of 
the target of interest, fragments generated from dynami¬ 
cal processes such as dissociation and Coulomb explosion, 
as well as background gases. The only way to separate 
out target specific electronic information is to correlate 
the detected electrons with their ionic partners [35J 03. 
An example of an ionic TOF detected after the ionization 
of Xe by intense 3.1 pm radiation is presented in Fig. [6] 
The singly ionized stable and long lived isotopes of Xe 
are resolved and their relative abundances agree with the 
known values |50| . Using a two particle coincidence con¬ 
dition, the electrons corresponding to the Xe + ions can 
be isolated from the rest of the electron background. In 
Sec. VIA we perform such an analysis to investigate the 
low momentum part of the electron spectrum and we 
present why many of the features evaded detection with 
TOF spectrometers. 


Statistics permitting, higher order particle coinci¬ 
dences are possible for multiple ionized ions. For ex¬ 
ample, the two peaks m corresponding to the doubly 
ionized 129 Xe isotope are visible in the inset of Fig. [6] As 


this is a double ion, it has two corresponding electrons. 
Using a three particle coincidence condition, we can iso¬ 
late the two electrons corresponding to each double ion. 
In Sec. |VI A| we present the results of such analysis and 
show, for the first time, an electron-electron correlation 
map resulting from NSDI triggered by intense rnid-IR 
radiation. 

The importance of coincidence detection becomes even 
more vital when molecular targets are investigated as 
they can undergo fragmentation upon SFI. In Fig. & 
the ionic TOF observed after the ionization of C 2 H 2 us¬ 
ing similar laser parameters is presented. Clearly many 
fragmentation channels can be observed. Without coin¬ 
cidence detection, electrons not associated with the frag¬ 
ment of interest become an unwanted background signal 
that can detrimentally affect data analysis or cloud inter¬ 
pretation. This point is discussed in detail with respect 
to performing LIED on polyatomic molecular targets in 
Ref. [37|. As SFP continues to progress from atomic 
and simple homonuclear diatomics towards more com¬ 
plex molecular targets, coincidence detection will become 
increasingly vital to extract specific information within a 
myriad of competing processes. 

There are currently two experimental methods that 
can achieve 3D detection in full coincidence: Firstly, ve¬ 
locity map imaging (VMI) |q 21 - |5T| is a technique where 
particles with the same mass and initial velocity are 
mapped onto the same position on a two dimensional 
(2D) detector. A VMI can provide electronic and ionic 
kinetic energies and angular distributions with full 47r sr 
collection whilst achieving energy resolutions on the per¬ 
cent level m- Most VMI apparatuses require cylindri¬ 
cal symmetry to transform the detected projection into 
a slice through the centre of the 3D momentum distri¬ 
bution, however, in special implementations this is not 
necessarily the case [56U591 . Moreover, VMIs typically 
only operate in full particle coincidence in special circum¬ 
stances m rather than by default. The second option is 
to use a ReMi [STHBB] , as is this case in this manuscript, 
which is based on cold-target recoil ion momentum spec¬ 
troscopy (COLTRIMS). The arrival time and position of 
all charged particles is detected in order to reconstruct 
the full 3D doubly differential cross-section of the inter¬ 
action. The operating principle of a ReMi is described in 
Sec. IVl A. 


V. EXPERIMENTAL SYSTEM 

Our experimental system consists of a high repetition 
rate mid-IR OPCPA source with long term power and 
carrier-envelope phase (CEP) stability in combination 
with a ReMi detection system. The mid-IR source pro¬ 
vides conditions in the QS regime while the ReMi images 
the entire 3D momentum distribution of the interaction 
in full coincidence. Overall, the system is currently ca¬ 
pable of creating and detecting electron energies over a 
six order of magnitude range and it provides the perfect 
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FIG. 6 : a) The measured ion TOF spectrum after the 
ionization of Xe using the 3.1 pm output. The 
individual Xe isotopes can be accurately resolved when 
low \E\ and \B\ extraction fields are used. The inset 
shows TOF information corresponding to 129 Xe 2+ . b) 
The same as in (a) except for C 2 H 2 and with high 
extraction fields. Many molecular fragments are 
observed with each one having its own corresponding 
electrons. 


conditions to investigate and utilize mid-IR SFI in the 
QS regime. 


A. Reaction microscope (ReMi) 

The functional principle of the ReMi detection sys¬ 
tem is as follows. The target gas is supersonically ex¬ 
panded into the ultra-high vacuum ( 10 ~ n mbar without 
gas load) and is subsequently skimmed in two succes¬ 
sive stages before reaching the ReMi. The expansion of 
the molecular beam thus features a decreased thermal 
motion of the particles and internal temperatures to the 
milikelvin range are feasible. Upon SFI, a combination 
of static electric (E) and magnetic (B) extraction fields 
direct all charged fragments to two opposing large-scale 
and position sensitive micro channel plate (MCP) detec¬ 
tors. The homogeneous electric field is created by a stack 
of metal electrodes inside the ReMi while the magnetic 


field is created by high current external coils. In Fig. 
the negative and positive particles created during the ion¬ 
ization event are represented by green and purple arrows, 
respectively. The 3D momentum vectors (with the x-axis 
along the laser propagation direction, the y-axis being 
anti-parallel to the direction of the gas jet and the 2 -axis 
along the direction of the static fields) and hence the ki¬ 
netic energy of the particles can be calculated from the 
measured TOF and the position of impact. The choice 
of the electric and magnetic fields define both the mo¬ 
mentum range that can be examined in the experiment 
and the achievable momentum resolution. By choosing 
relatively low fields, high momentum resolutions can be 
achieved, as presented in F ig. [T)r. Here the total mo¬ 
mentum resolution (A p = y' Ap* + A p 9 + A pi) has been 
calculated [S3] as a function of the longitudinal momen¬ 
tum (p|| = p z ) while \E\ = 1.3 V/cm is kept constant 
for both \B\ = 4.4 G (red curve) and \B\ = 5.1 G (blue 
curve). The observed discontinuities for both magnetic 
fields are due to electrons with integer multiples of the 
cyclotron frequency hitting the same position on the de¬ 
tector. While these so-called ‘nodes’ seem inconvenient at 
first, they turn out to be useful as they serve as an accu¬ 
rate and straight forward measurement of B. Moreover, 
is it easy to remove these regions by taking two measure¬ 
ments with different magnetic fields and combining the 
areas with the highest momentum resolution. The results 
of this technique are represented by the grey shaded re¬ 
gion in Fig. 03 where A p < 5 x 10 ~ 2 a.u. for all values of 
p\\. In terms of energy these values correspond to detec¬ 
tion resolutions of several meV. It is straightforward to 
further reduce the electric and magnetic fields to achieve 
even better momentum resolution. By choosing higher 
values of \E\ and \B\ it is possible to detect re-scattered 
electrons with kinetic energies of up to 1 keV. The magni¬ 
tude of E controls the maximum detectable longitudinal 
electron momentum component (jp\ |) while the maximum 

transverse component (pj_ = \ pi + p 2 ) is dependent on 

the magnitude of B. Figure 0: presents a logarithmi¬ 
cally scaled 3D momentum distribution detected under 
these high static field conditions (\E\ = 51 V/cm and 
\B\ = 39 G) where electron energies of > 500 eV were 
detected. 

B. Mid-IR OPCPA system 

The OPCPA system [BH |B5] is based on a multi¬ 
color fiber frontend with difference frequency generation 
(DFG) between the two outputs supplying broadband 
mid-IR pulses at 3.1 ;im center wavelength. The DFG 
stage is a crucial component of the OPCPA as it provides 
passive carrier-envelope phase (CEP) stability optically. 
The intrinsic CEP stabilization is a decisive advantage 
over electronically stabilized systems which manifests it¬ 
self in unsurpassed long time stability of our system. 











(a) 



FIG. 7: a) The experimental apparatus used to study SFP at long wavelengths. The operating principle of the ReMi 
is described in the text. The 160 kHz OPCPA has two outputs that can be focussed into a supersonically expanding 
gas jet to induce SFI. The full 3D momentum distribution of the interaction can be mesaured. b) The momentum 
resolution of the ReMi as a function of p\\ for two low values of the magnetic field. The ‘nodes’ can be eliminated by 
stitching the results of multiple measurements with different values of B together. 


The system provides both 1.7 /im and 3.1 ^m radiation. 
At 3.1 /im, pulse durations below three optical cycles 
FWHM are achieved (32 fs = 2.9 optical cycles) [66]. The 
3.1 /im output is passively CEP stable and a stability of 
250 rnrad rms over 11 minutes has been measured |64| . 
Without any additional effort, pulse durations of 100 fs 
are reached for the signal at 1.7 /im. The 15 /iJ and 
20 /iJ pulse energies of the 1.7 /im and 3.1 /im wave¬ 
lengths, respectively, are focused with an on-axis, gold- 
coated paraboloid (focal length of 50 mm) into the gas 
jet of the ReMi. Intensities in the range of 10 14 W cm -2 
are readily achieved for both wavelengths. 


C. The importance of repetition rate and stability 

A major challenge in performing SFP experiments with 
long wavelength sources is quantum diffusion m ■ Here 
the electron wavepacket experiences larger longitudinal 
and transverse spreading due to the longer excursion time 
at increased wavelength. Consequently, the larger extent 
of the recolliding wavepacket results in a smaller recom¬ 
bination or recollision cross-section which manifests in 
reduced count rate when driving SFP at longer wave¬ 
length mm- Experimental evidence shows that the 
ratio of re-scattered to directly emitted electrons scales 
roughly as A -4 [8|. 

We choose to highlight this crucial point for experi¬ 
mental investigations in Fig. [8] by defining an arbitrary 
standard unit of ‘acquisition time’ (AT) to be the time 
required to reach a certain signal-to-noise ratio of re¬ 
scattered electrons when using a standard 1 kHz, 0.8 /im 


radiation source (dashed black line). Using the above 
mentioned A” 4 scaling and assuming all other experimen¬ 
tal parameters stay constant, we can estimate the number 
of standard ATs required to reach the identical signal-to- 
noise ratio as a function of driving wavelength. The re¬ 
sults of this calculation are presented in Fig. [8] for 1 kHz 
sources (blue curve). The triangular data points repre¬ 
sent the wavelength of some sources that have recently 
been used to study electron re-scattering mmm- The 
required ATs for identical signal-to-noise scale rapidly 
with wavelength; compared with 0.8 /mi radiation we 
predict a one and two orders of magnitude increase for 
1.4 /im and 2.5 /im radiation, respectively. At a wave¬ 
length of 3.6 /im the required AT is over 400. We would 
like to put the arbitrarily defined ATs and the quoted 
numbers into perspective: if a standard AT is equal to 
one hour then at 3.6 /im the new experimental time would 
be over 17 days. Clearly, the unfavorable wavelength 
scaling of SFP must be addressed to ensure continual 
advancement of the field. 

An obvious solution to this dilemma is to increase 
the repetition rate of the driving laser [HI I68U70] , 
Ti:Sapphire based laser systems u n m sa are typi¬ 
cally limited to a few kHz; the range of acquisition times 
that can realistically be achieved with such systems is 
represented by the blue shaded region in Fig. [8] Mid- 
IR sources at much higher repetition rates, i.e. above 
several tens of kilohertz, can be realized with a range of 
techniques but so far only OPCPA has provided a viable 
solution. High average power fiber laser sources are oper¬ 
ational, but additional down conversion into the mid-IR 
carries a dramatic penalty in pulse energy. Moreover, 
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FIG. 8 : The number of ‘acquisition times’ (as defined in 
the text) required for a range of driving wavelengths 
after accounting for the A -4 scaling of electron 
re-scattering. Typical 1 kHz systems based on 
Ti:Sappliire technology (Ref. [IT) - blue triangles) are 
limited as total acquisition times become untenable. 
Even the expected repetition rate upper limit of 5 kHz 
isn’t fast enough to compensate for the decrease of 
signal (shaded region). Using a 160 kHz system such as 
our OPCPA, on the other hand, allows much shorter 
acquisition times (purple circles). 


few-cycle durations and CEP stability are not possible 
with such sources at high intensities (> 10 13 Wcm -2 ). 
Another approach could be based on passive enhance¬ 
ment cavities m in the mid-IR but such apparatuses are 
extremely challenging and are yet to be implemented out¬ 
side of the near-IR. Due to these technical problems we 
have devised a radically new platform mm (summa¬ 
rized in Sec. V B) that leverages the concept of OPCPA to 
achieve the combination of mid-IR wavelength (3.1 /im), 
few-cycle pulse duration (> 2.9 cycles) with CEP stabil¬ 
ity, SFP relevant intensities (< 3 x 10 14 W cm -1 ), 160 
kHz repetition rate and long term stability (< 1% rms 
over 4.5 hours) [55J. Additionally, parasitically generated 
wavelengths such as the OPCPA’s signal at 1.7 /irn can 
be used for further frequency down-converison. When 
the above ‘acquisition time’ analysis is applied to the two 
wavelengths of our system it is found that the required 
experimental time is actually an order of magnitude less 
than a standard AT for 1.7 /irn and only 40% longer for 
3.1 /iin (the purple circles in Fig. [ 8 ]). Combining these 
points with the fact that the mid-IR source can actually 
reach a 7 < 0.2, it is clear that mid-IR OPCPA solves 
the problem of studying SFP at longer wavelengths. 


VI. EXPERIMENTAL CAPABILITIES 

We now demonstrate the capabilities of our approach 
by investigating strong-field electron recollision at low 
and high energies. We choose Xe as an example for which 


we measure a 550 eV wide kinetic energy spectrum after 
ionization with intense 3.1 fim radiation. By optimising 
the electric and magnetic guiding fields of the ReMi, as 
discussed in Sec. [VJ we investigate different parts of the 
spectrum. With low fields we resolve meV features at 
energies below 1 eV while with high fields we collect all 
ejected electrons up to energies of hundreds of eV. 


A. Low energy recollisions 


As already described briefly in Sec. [TJ experimental in¬ 
vestigation of SFI with longer wavelength laser sources 
(A > 1.5 /im) revealed unexpected behavior in the ob¬ 
served electron energy spectra that deviated from SFA 
predictions urn- In addition to the well-known form of 
the kinetic energy spectrum and the drop in yield beyond 
directly emitted electrons, a peaked structure at low en¬ 
ergies (1 — 3 eV) was observed. This feature, dubbed the 
low-energy structure (LES), was found to result from the 
interplay between the Coulomb binding potential and the 
laser electric field [241 - B6] . Additional features have been 
observed at even lower energies and were dubbed very low 
energy structures (VLES) [23j. I n these cases the descrip¬ 
tion of the interaction using SFA based models seemed 
to be missing crucial aspects of the post ionization be¬ 
havior of the electron in the laser field. It was found 
that multiple ‘revisits’ of the electron close to the ionic 
core, and therefore within the influence of the Coulomb 
field, had the affect of bunching the oscillating electrons 
in momentum space l‘2?]j.‘{2| . It is interesting to note that 
classical trajectory calculations that ignore the Coulomb 
potential were able to reproduce the LES features [Z2j, 
however, the Coulomb potential markedly enhances the 
effect. From an experimental point of view it is important 
to remark that the LES was originally measured using a 
TOF spectrometer, which, as explained in Sec. 0 only 
detects a small fraction of the entire momentum distribu¬ 
tion. This restriction in detection translated into incom¬ 
plete coverage of the momentum space, thereby missing a 
large fraction of information from the observed effects. A 
comprehensive picture of the observed effects was gained 
from our full 3D and high resolution momentum measure¬ 
ments with a ReMi mm- In these measurements, an 
additional feature near zero momentum was observed and 
called the zero energy structure (ZES). Recently, we per¬ 
formed a combined experimental and theoretical investi¬ 
gation }35| in which the origins of all three features were 
quantitatively confirmed to be due to two-dimensional 
Coulomb focusing (for the LES and VLES) and post¬ 
ionization recapture of electrons into high-lying Rydberg 
states with subsequent field ionization (for the ZES). The 
latter mechanism is also known as frustrated tunnel ion¬ 
ization mm- 
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FIG. 9: a) The measured electron momentum 
distribution corresponding to Xe + after interacting with 
an intense 71 fs, 3.1 pm pulse with a peak intensity of 
I = 4 x 10 13 W cm -2 (7 = 0.4). The data has been 
symmetrized with respect to py = 0. Line-outs at three 
values of p± are presented in (b) to highlight the 
different LES orders and the VLES (black dashed 
circles). The white dashed line in (a) represents the 
acceptance cone of a TOF spectrometer with ©a = 5°, 
which roughly precludes the detection of the features 
within the grey shaded region in (b). 


1. Singly ionized Xe 


Figure [9] presents kinematically complete measure¬ 
ments of single ionization of Xe using 3.1 pm radiation at 
a peak intensity of I = 4 x 10 13 W cm -2 and a pulse du¬ 
ration of 71 fs (7 optical cycles, 7 = 0.4). In Fig. the 
resulting electron momentum distribution corresponding 
to Xe + is presented in cylindrical coordinates after inte¬ 
gration over the azimuthal angle (see Fig. 11 for further 
details). While the symmetric distribution shows similar 
characteristics to studies performed on other targets such 
as Ar [H [Ml 133, o 2 HU and N 2 [M], more low energy 
features seem to be visible. We find a broad distribution 
of direct electrons with pronounced low energy features, 
with the most distinct being the merged VLES and ZES 
near p|| ss 0 a.u. and p± < 0.03 a.u.. Note that com¬ 
pared to Ar, O 2 and N 2 , the ZES in Xe does not appear 
as clearly (see Fig. 2 of Ref. [M] for a direct compari¬ 
son). The experimental conditions for Xe, compared to 
the other species we investigated, included a lower peak 


intensity of the laser (/ = 4 x 10 13 W cm 2 measuring 
Xe compared to I = 9 x 10 13 W cm -2 measuring Ar ESI) 
despite a higher extraction field of the ReMi (2.0 V/cm 
in the case of Xe compared to 1.3 V/cm in the case of 
Ar ESI)- To highlight the individual features three nor¬ 
malized horizontal line-outs are presented in Fig. HJo at 
values of p± = 0.02 a.u. (red curve), p± = 0.03 a.u. 
(blue curve) and p± = 0.07 a.u. (orange curve). The 
p_ l = 0.03 a.u. and p± = 0.07 a.u. curves have been 
offset by +0.2 and +0.4, respectively, to enhance visibil¬ 
ity. The p l = 0.03 a.u. curve shows the richest structure 
with five different structures being resolvable: LESi at 
p\\ ~ 0.60 — 0.70 a.u., LES 2 at py ss 0.40 a.u., LES 3 
at py ~ 0.25 a.u., LES 4 at py « 0.15 a.u. The charac¬ 
teristic V-shape associated with the VLES [53] is visible 
between |py | < 0.06 a.u. The positions of the LES orders 
(dashed black circles) agree well with predictions [291 . 
For the p_L = 0.02 a.u. and p_L = 0.07 a.u. line-outs 
the visibility of the various LES orders vary. The other 
most striking feature is the centre peak (py = 0) for 
pj_ = 0.02 a.u., which corresponds to the ZES. Note that 
the white dashed line in Fig. §1 represents the acceptance 
cone for a TOF spectrometer with 9 a = 5° and with its 
axis along the polarization direction. Higher transverse 
momentum parts of the LES orders, the VLES and the 
ZES, such as those within the grey shaded area in Fig. n+ 
are not detected when using such an apparatus. This elu¬ 
cidates the need for 3D imaging of the momentum space 
in order to observe all low energy recollision features. 


2. Doubly ionized Xe 


Next, we use the coincidence detection capabilities to 
extract electrons correspo ndi ng to doubly ionized xenon 
atoms Xe 2+ only. Figure 10 shows the first coincidence 
measurements of doubly ionized Xe with mid-IR fields, 
and we obtain high measurement statistics due to the 
stability and repetition rate of our OPCPA system. The 
data exhibits a distribution similar to the Xe + data with 
three possible structures seemingly emerging above the 
noise floor: 1 ) a feature at py = ± 0.6 a.u. and pj_ « 
0.07 a.u., which is close to LESi for Xe + ; 2) a feature near 
py = ± 0.3 a.u. andpj_ w 0.04 a.u., which is close to LES 2 
or LES 3 for Xe + ; and 3) a VLES like distribution at py = 
0 a.u. and pj_ ~ 0.02 a.u.. These results indicate, for 
the first time, that electrons also could undergo multiple 
recollisions with the ionic Coulomb field during double 
ionization. 

To determine the mechanism behind mid-IR double 
ionization [73 EH we plot the parallel momentum distri¬ 
bution of the double ions (see Fig. llOb) and the cor¬ 
relation between the two electrons Re) ,e 2 ). Fig. 10:). 
The parallel momentum distribution of the double ion, 
in Fig. clearly shows the double hump structure 

which is characteristic of NSDI EZI and is not explained 
by sequential double ionization (SDI), where two elec¬ 
trons are emitted independently at subsequent laser field 
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tions parallel to the horizontal and vertical axes of the 
(ei , 62 ) correlation plot. Figure 10; shows features remi¬ 
niscent of both (e,2e) and RESI. Interestingly, upon first 
inspection, the shape of the measured correlation shows 
a striking similarity to Fig. 7c of Ref. [3U] where RESI 
was found to be the dominant mechanism for very dif¬ 
ferent experimental parameters. In that paper, Bergues 
et al. [80 observed similar features in Ar for near-single¬ 
cycle pulses at A = 750 nm and at I = 3 x 10 14 Wcm” 2 . 
Our conditions for Xe are markedly different with a 6.5 
cycle pulse (71 fs) at I = 4 x 10 13 W cm ” 2 in the mid-IR. 
Without accompanying theoretical calculations it is diffi¬ 
cult to determine the exact contribution of either mech¬ 
anism, however, theoretical investigations into mid-IR 
double ionization are only now starting to become avail¬ 
able [5TH53] . When combined with these first ever full 
coincidence mid-IR double ionisation results theoretical 
models will be able to unravel the mechanisms of double 
ionisation within the mid-IR regime. 


B. High energy recollisions 


FIG. 10: a) The electron momentum distribution 
corresponding to Xe 2+ ions for the same experimental 
conditions as in Fig. b) The longitudinal momentum 
distribution of the X? + ions (blue curve) where a dip 
at p |i = 0 indicates the involvement of electron 
recollision. For reference the single peak distribution of 
the Xe + ions is also presented (black curve) along with 
a Gaussian fit (grey curve), c) The correlated 
two-electron momentum distribution, which features a 
cross-like shape with a preference to the first and third 
quadrants. 


maxima. NSDI exhibits two mechanisms: 1) the recol¬ 
liding electron can cause the direct emission of a second 
electron via electron-impact ionization, otherwise known 
as (e,2e) ionization [19] ; or 2) recollision excitation with 
subsequent ionization (RESI) [73] where the first elec¬ 
tron excites the parent ion which leads to tunnel ion¬ 
ization of a second electron. The form of the double 
hump in Fig. | 10 [ 3 , i.e. the central dip not being deeply 
pronounced, is already an indication that both mecha¬ 
nisms of NSDI are present. Further information can be 
extracted from the electron correlation plot in Fig. [To] c. 
The spectrum is symmetric about P|| el = P ||, e 2 due to 
the indistinguishability of the two electrons in the exper¬ 
iment. To interpret the measured correlation pattern, 
we recall that a correlation (first and third quadrants) 
means that both electrons leave in the same direction 
(and in the direction opposite to the double ion) f79] . 
while an anti-correlation (second and fourth quadrants) 
means that the electrons escape in opposite directions to 
each other (the double ion exhibits no momentum) [76] . 
The RESI mechanism was described previously [ 73] and 
due to energy sharing exhibits ’cross-shaped’ distribu- 


The capability of our apparatus towards investigat¬ 
ing SFI at high electron energies is demonstrated next. 
We measure doubly differential electron scattering cross- 
sections (DCSs) using LIED [51H53] . The realization of 
LIED is dependent on two key points: 1) the genera¬ 
tion of high impact energies for the recolliding electrons; 
and 2) the ability to extract field free DCSs despite the 
presence of the strong laser field. Both conditions are 
fulfilled with our apparatus. Previously, such measure¬ 
ments of the DCS were either limited due to the shorter 
wavelength of the Ti:Sapphire laser, or the inability of 
collecting the full 3D momentum distribution mi (571- 
eh. Here, we show the first LIED experiment based on 
3D momentum detection in full particle coincidence. 

Mid-IR LIED is achieved by operating the ReMi at 
high static electric and magnetic fields to detect the 
entire 3D momentum distribution of the emitted elec¬ 
trons. Figure |ll[ t displays the momentum distribution 
created during the interaction of our intense 3.1 pm ra¬ 
diation with Xe atoms. After an integration over the 
azimuthal angle (f) and the application of the appropriate 
co-ordinate transformation (Jacobian J = l/|pjJ), a 2D 
momentum map representation of the detected momen¬ 
tum distribution is extracted (Fig. m- The momentum 
map covers a six order of magnitude range in counts and 
features an unprecedentedly rich structure. DCSs are ex¬ 
tracted from the 2D momentum distribution (or diffrac¬ 
tion pattern) in a manner similar to conventional elec¬ 
tron diffraction (CED) provided that three critical crite¬ 
ria are satisfied: 1) QS conditions must be achieved so 
that the final momentum distribution of electrons can be 
expressed as the product of the field-free DCS multiplied 
by the returning electron wave packet [92]. In this case, 
the vector potential of the driving field can easily be cal¬ 
culated and taken into account and field-free DCSs can be 
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FIG. 11: The extraction of elastic DCSs of electrons scattering off Xe atoms using the LIED technique, a) Xe is 
ionized by the 3.1 //in radiation with an intensity of~ 6x 10 13 W cm -2 and the full 3D momentum distribution is 
detected by the RcMi. b) After integration around the azimuthal angle (cf>) and the application of the appropriate 
co-ordinate transformation Jacobian a 2D representation of the momentum distribution is achieved, c) The DCSs 
are found by tracing the radius (corresponding to the return electron momentum fc r ) of a circle that has been shifted 
by the vector potential A r at the time of electron rescattering. Accurate DCSs are presented for return energies E T 
of 75 eV (blue circles), 100 eV (red squares) and 125 eV (green triangles) over a wide scattering angle range. The 

solid lines are taken from the NIST database [55] , 


extracted; 2) the returning electron wavepacket must be 
energetic enough to probe the core of the target and not 
the valence electronic structure (i.e. E r > 50 eV); and 
3) the recolliding electron must propagate sufficiently far 
away from the interaction region such that it approxi¬ 
mates a plane wave upon return. These three require¬ 
ments can only be simultaneously satisfied by intense 
mid-IR sources (I > 5 x 10 13 W cm -2 and A > 2 //m). 
Once these stringent conditions are satisfied, the achieved 
momentum transfers are similar between LIED and CED 
despite the difference in incident electron energies and 
scattering angles. That means in mid-IR LIED, elec¬ 
tron energies E > 50 eV that scatter at angles between 


30° < 0 r < 180° are typically utilized while in CED the 
electron energies are on the order of 10 4 —10 5 eV and scat¬ 
ter at forward angles (0° < 0 T < 10°). This means that 
the classical momentum transfer (q = 2fc r sin$ r /2, where 
fc r is the electron incident momentum) is comparable for 
both methods and, as such, mid-IR strong-field induced 
re-scattered electrons contain structural information of 
the target gas. 

The theoretical framework behind the analysis of the 
mid-IR LIED pattern is provided by the quantitative re¬ 
scattering (QRS) theory I92U94I . This model states that 
angularly resolved DCSs (up to an overall scaling factor) 
can be extracted from the measured electron momentum 
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distribution by tracing the measured counts along the cir¬ 
cumference of a circle of radius fc r that has been shifted 
from the origin by the vector potential of the laser at the 
time of re-scattering (A r ). Here fc r = y/2E r is the inci¬ 
dent momentum of the re-scattering electron with energy 
E t . The relation between the incident energy of the re¬ 
turning electron E T and the final electron energy of the 
rescattered electron E^ n is plotted in the inset of Fig. [l] 
for both long and short trajectories (derivation from e.g. 
in section II.A. of [92]). An example of this extraction 
procedure for 75 eV (blue), 100 eV (red) and 125 eV 
(green) energy electrons is presented in Fig. [H}> The 
extracted DCSs are presented in Fig. ES as a function 
of the scattering angle. As each returning electron en¬ 
counters a slightly different vector potential, the rescat¬ 
tering circles are shifted by differing amounts. The 75 eV 
(blue circles) and 125 eV (green triangles) data have been 
scaled by factors of 10 1 and 10“ 1 compared to the 100 
eV (red circles) data, respectively, in order to enhance vi¬ 
sualization. The error bars are derived from Poissonian 
statistics and are propagated through the integration and 
co-ordinate transformation procedures. The large error 
bars at 9 r ss 180° are due to the application of the Ja¬ 
cobian mentioned above. Also presented in Fig. Ej} are 
the Xe DCSs provided by NIST [95], which have been 
scaled to overlap with the experimental data. Our data 
are able to accurately resolve most of the predicted rich 
structure for 9 r > 60°. For these angles the maxima and 
minima as well as slope changes are generally well re¬ 
produced while for smaller scattering angles ( 0 r < 60°) 
the data is ’contaminated’ through background arising 
from direct electrons. While these direct electrons cause 
a slight mismatch between the experimental and theoret¬ 
ical data, the minima between 9 r = 50 — 60° are still re¬ 
solvable. Overall we observe an excellent match between 
our measured data and the theoretical expectations. We 
note that the presented NIST data is calculated for neu¬ 
tral Xe. It has been shown that the electron energies uti¬ 
lized here (FAm > 50 eV) are core penetrating and can 
be used to extract structural information m ■ Therefore 
the electronic structure of the target has a negligible ef¬ 
fect on the extracted DCSs for these experimental condi¬ 
tions. These results show that field free DCSs can indeed 
be extracted over a wide range of energies and scattering 
angles at unprecedented accuracies. 

VII. SUMMARY AND OUTLOOK 

We have presented a unique methodology to investi¬ 
gate mid-IR driven SFP with a high repetition rate mid- 
IR OPCPA system in combination with a ReMi that is 
able to detect ionization fragments over a six order of 
magnitude energy range. This apparatus now permits 


unambiguous investigations of SFI and related effects in 
the deep QS regime, thereby allowing experimental obser¬ 
vations to be scrutinized against theoretical descriptions. 
We demonstrate our concept across the entire accessi¬ 
ble kinetic energy scale by investigating SFI and electron 
recollision with Xe atoms. We observe low energy struc¬ 
tures up to the fourth order, the VLES and the ZES 
in Xe but find a much reduced ZES that would arise 
from capture into Rydberg states and subsequent ReMi 
field ionization. This surprising absence of a strong ZES 
peak in Xe already warrants further investigation as it 
will prompt information about the FTI process. Next, 
first coincidence measurements of mid-IR driven NSDI 
are shown. The analysis of the data hints at the com¬ 
bined action of (e,2e) and RESI mechanisms. Interest¬ 
ingly, this very first (el,e 2 ) correlation map in the deep 
QS regime (7 = 0.4) for a 6.5 cycle pulse exhibits fea¬ 
tures that were observed for entirely different conditions 
with a near-single-cycle CEP-stable pulses at A = 750 nnr 
in Ar and for 7 = 0.7 [5U|. Lastly, we exploit the high en¬ 
ergy detection capabilities of our apparatus with mid-IR 
LIED. We achieve core penetrating conditions that per¬ 
mit extraction of the doubly differential elastic electron 
scattering cross section in Xe and the retrieved DCSs 
match published NIST data very well. Probing atomic 
structure with mid-IR LIED carries clear potential for 
the realization of femtosecond time resolved imaging of 
molecular structure and for tracking atomic constituents 
during structural changes such as isomerization or disso¬ 
ciation. In fact, recent results from our group show that 
the method can indeed be used to determine the entire 
conformation of polyatomic acetylene (C 2 H 2 ) [37]. 
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